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ABSTRACT 

Until today, the small size of magnetic elements in quiet Sun areas has required the application of 
indirect methods, such as the line-ratio technique or multi-component inversions, to infer their physical 
properties. A consistent match to the observed Stokes profiles could only be obtained by introducing a 
magnetic filling factor that specifies the fraction of the observed pixel filled with magnetic field. Here, 
we investigate the properties of a small magnetic patch in the quiet Sun observed with the IMaX 
magnetograph on board the balloon-borne telescope Sunrise with unprecedented spatial resolution 
and low instrumental stray light. We apply an inversion technique based on the numerical solution of 
the radiative transfer equation to retrieve the temperature stratification and the field strength in the 
magnetic patch. The observations can be well reproduced with a one-component, fully magnetized 
atmosphere with a field strength exceeding 1 kG and a significantly enhanced temperature in the mid- 
to upper photosphere with respect to its surroundings, consistent with semi-empirical flux tube models 
for plage regions. We therefore conclude that, within the framework of a simple atmospheric model, 
the IMaX measurements resolve the observed quiet- Sun flux tube. 

Subject headings: Sun: magnetic topology — Sun: photosphere — techniques: polarimetric — tech- 
niques: spectroscopic 



1. introduction 



Ever since the 1970 s (e.g. iHoward fc Stenflol [1972; 
Fr azier fc Stenflo 1972|) it has been known that small- 
scale magnetic elements are unresolved. Therefore, the 
magnetic field strength averaged over the resolution el- 
ement was found to be significantly lower than the field 
strength of the magnetic structure. Their intrinsic prop- 
erties (field strength and temperature structure) could 
only be obtained by introducing a magnetic filling factor 
to account for the fact that only a certain fraction of the 
observed resolution element contains the magnetic struc- 
ture. This factor is sometimes also called a stray light 
factor, although the physi cal meaning is somewhat differ- 
ent. Indirect methods (see Solanki 1993, for an overview) 
such as the line-ratio techni que ([St enflo 19 731) or the in- 
version of line profiles (e.g. Martinez Pil let et "all Il997t ) 
pointed to the kilo-Gauss nature of the small-scale mag- 
netic elements in the network. The reality of these kilo- 
Gauss field strengths has been confirmed by direct split- 
ting of infrar ed lines (iHaryev Hall ll975l: lHarveylll977l 
iRabinl 119921 : iRiiedi et all 11992ft. Bv employing speckle 
imaging techniques, Keller (1992) observed Stokes V am- 
plitudes in plage regions close to the values predicted 
by semi-empirical flux tube models. This groundbreak- 
ing analysis suffers from the assumption that the chosen 
empirical atmosphere is applicable to the particular ob- 
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served feature. 

Another basic property of small magnetic elements 
is that they are hotter than their surroundings in the 
middle and upper photosphere, with intrinsic tempera- 
ture difference s of ^1000 K being reached at logr ~ —2 
(Solanki 1986). These high temperatures lead to a weak- 
ening of many Fraunhofer lines, resulting in lin e gaps or 
discontinuities on spectrograms (jSheelevlll967[ ). 

Here, we report on the self-consistent determination of 
the magnetic field strength and temperature in the mag- 
netic elements, without taking recourse to a magnetic 
filling factor. Whether the observed magnetic elements 
have finally been resolved can then be judged by com- 
paring with both the field strengths and the excess tem- 
perature in magnetic elements given in the literature as 
deduced with methods that allow for the magnetic filling 
factor. 

2. OBSERVATIONS 

We use dat a from the Imaging Magnet ograph exper- 
iment (IMaX, iMartmez Pillet et all 1201 Of) on b oard the 
Sunrise balloon mission ([Barthol et al.M201Q[ ). For an 
overview of the Sunrise d ata set and a descri ption of 
selected results we refer to [Solanki et al.l (|2010h . IMaX 
measured the full Stokes vector in five wavelength po- 
sitions located at (-80, -40, +40, +80, and +227) mA 
from the line center of the magnetically sensitive Fel 
5250.2 A line (Lande factor g = 3). The data set was 
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Fig. 1. — IMaX continuum intensity map (top) and Stokes V 
map (bottom) normalized to the continuum level of the average / 
profile (see Figure [2]). The box encloses the region for the inversion 
maps in Section 0] 

obtained on 2009 June 9, at 00:36:45 UT (data set 163- 
209). The observed region was located close to the disk 
center. The total field of view of IMaX covered 50"x50" 
and contained a quiet- Sun area with isolated network el- 
ements, visible in the Stokes V map. The Stokes V map 
was obtained by integrating over the two blue wavelength 
points (—80 and —40mA from the line center). For the 
analysis in this Letter we discuss in detail a 7 ."7 x 4 ."7 
large subfield containing two small network patches with 
a dimension of «1" (see Figure [TJ position x =10. // 5, 
y =5. "8 and x =12. "5, y =8. // 3). On average we find 
roughly 10 network patches with similar Stokes signals 
per IMaX snapshot. 

The standard IMaX data reduction routines were used 
to perform dark-current subtraction, flat-field correc- 
tion, and cross-talk removal. The noise level in the 
Stokes V, Q and U signals was ^ 3 x 10 -3 . A spa- 
tial resolution of Of'15— 0."18 in all Stokes parameters 
was achieve d by combining the on-board image stabiliza- 
tion system (jGandorfer et alJ201ttlBerkefeld et al.ll2010D 
with a phase-diversity -based post processing technique 
([Martinez Pillet et all [201 0( 1. 

3. SETUP OF INVERSIONS 

The inversions of the Stokes vector observed with 
IMaX are carried out using t he SPINOR code 
(jFrutker et al.l 120001 : iFrutkerl l2000l ) . This code numer- 
ically solves the radiative transfer equation (RTE) un- 
der the assumption of local thermodynamic equilibrium 
and minimizes the difference between the measured pro- 
file and the computed synthetic profile using a response- 
function-based Levenberg-Marquardt algorithm. The 
limited number of wavelength points of the IMaX mea- 
surements requires the usage of a simple model at- 
mosphere with height-independent values for the mag- 
netic field vector and the line-of-sight velocity. This 



model atmosphere is in some ways comparable to Milne- 
Eddington atmospheres that involve nine free parame- 
ters, four of them being so-called ad hoc parameters 
without a one-to-one connection with atmospheric pa- 
rameters (i.e., the line center to continuum opacity ra- 
tio, two parameters describing the linear source function, 
and damping). In contrast, our numerical solution of the 
RTE requires only eight free parameters: the magnetic 
field strength, £>, inclination, 7, azimuth, x-> the line-of- 
sight velocity, ^los, macro- and micro-turbulence, £ mac 
and £mic> and two parameters describing the temperature 
stratification. We used the temperature str atification 
defin ed by the HSRASP model as the basis (| Chapman! 
|1979[ ) and modified it by adding a linear function de- 
fined by an offset, To, and a temperature gradient offset 
(changing the gradient in the HSRASP profile), TQ r ad : 

T(logr) = T HS RASp(logr)+To + T Grad .logr. (1) 

Positive values of Tcrad lead to a steepening of the tem- 
perature profile with respect to HSRASP while negative 
values lead to a flatter profile. To avoid a temperature 
increase for large negative values of Tcrad in the upper 
atmospheric layers, we forced the gradient of T(logr) to 
be steeper than 150 K per unit of logr, similar to the 
gradie nt of the PLA atmosphere of Solanki & Bri glievicl 
(| 19921 see Section 3]) at a height of logr = — 3. As a con- 
sequence of the high spatial resolution of the IMaX obser- 
vations, we set the macro-turbulence velocity to 0. For a 
proper comparison with the measured IMaX profiles, we 
convolve the synthetic profiles with the filter curve, mod- 
eled with a Gaussian with a full widt h at half- maximum 
of 85 mA (M artinez Pillet et al.ll20Toh . 

In Figure El the I profile computed from the HSRA 
model is overplotted on the profile averaged over the 
IMaX field of view. This average quiet-Sun profile (dia- 
mond symbols in Figure [2]) is consistent with a synthetic 
profile (solid red line) computed from the HSRA model 
(jGingerich et al.l Il97l[ ) assuming a micro- and macro- 
turbulence of 0.8kms _1 and Okrns" 1 , respectively, and 
then convolved with a 85 mA Gaussian, corresponding 
to the bandpass of the IMaX etalon. The good match of 
the measured average Stokes I profile with the synthetic 
profile computed from the HSRA atmosphere justifies 
the usage of this atmospheric model as the basis for the 
inversions. The temperature of the HSRA atmosphere 
at logr = is 6390 K with a gradient at this height of 
1900 K per unit of logr. We normalize all Stokes profiles 
to the average quiet- Sun continuum level. This normal- 
ization replaces the photometric calibration necessary to 
use Planck's law for determining the temperatures in the 
continuum forming layers. The intensity in the line core 
constrains the temperature in higher layers of the solar 
atmosphere. For later comparison (see Section [4j), we 
overplot with cross symbols the I profile measured in 
one of the network patches marked with (1) in Figure [IJ 
normalized to the average quiet- Sun continuum. 

Neighboring spectral lines might be shifted to the 
IMaX wavelength interval when large line-of-sight ve- 
locities are present. In this analysis, we therefore in- 
cluded the Col line at 5250.0 A and the Fel line at 
5250.6 A These lines are calculated assuming the same 
atmospheric model as used for the synthesis of the IMaX 
line (Fel 5250.2 A). The decrease in intensity at both 
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Fig. 2. — Stokes I profiles. The diamonds represent the profile 
measured at the five IMaX wavelength points averaged over the 
field of view observed with IMaX. The solid line is a synthetic 
profile resulting from the HSRA model. The crosses correspond 
to the measured / profile in a network patch marked with (1) in 
Figure [T] 

boundaries in the / profiles shown in Figures [2] and |4] 
stems from the neighboring Col and Fel lines. 

4. RESULTS 

The inversions were applied to an IMaX set of Stokes 
parameter images. We now discuss the results obtained 
in a box of 7 "7 x 4. // 7 (140 x 85 pixels), containing a 
region of enhanced Stokes V signal (see Figured]). The 
resulting maps for the magnetic field strength (constant 
with height), the temperature at logr = +0.5 and —2.0, 
and the temperature gradient offset Tcrad (see Equa- 
tion [T]) are shown in Figure G3 Contour lines of the mag- 
netic field strength are overlaid on the temperature and 
temperature gradient maps, contour lines of the temper- 
ature at logr = — 1.0 are plotted on top of the magnetic 
field strength map. The magnetic field outside the mag- 
netic patches is weak and lies close to or below the detec- 
tion limit of IMaX. In the center of the magnetic patch 
(1) the field strength reaches values of up to 1.45 kG. 
Note that the inversions return these field strengths with- 
out introducing a magnetic filling factor. This is one 
among many other examples already found in the IMaX 
data set. 

The temperature map at logr = +0.5 (panel (b) of 
Figure [3]) clearly reflects the typical structure of the 
quiet- Sun in the deep layers of the photosphere, with 
hot granules and cool inter granular lanes (see, e.g., 
Bor rero fe Bellot Rubiol l200l . With increasing height, 
the temperature in the granules falls off more rapidly 
than in the intergra nular lanes, resulting in a reversed 
granulation pattern (jCheung et al.l 120071 ). Panel (c) of 
Figure [3] shows the temperature at logr = — 2.0, a height 
between these extremes. At this height, hotter regions 
generally have stronger magnetic fields, indicating a shal- 
low temperature gradient in the magnetic regions. This 
is reflected in the map of the temperature gradient off- 
set, Tcrad, shown in panel (d) of Figured In the gran- 
ules the temperature falls off more rapidly than in the 
HSRASP atmospheric model (blue colors), whereas in 
the intergranular lanes, and especially in the magnetic 
field patches, the gradient is significantly lower (negative 
values of Tcrad)- As a result of this low gradient, the 
temperature in the kilo- Gauss features at logr = — 2.0 
is ^1000 K higher than in the granules surrounding the 
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Fig. 3. — Maps of magnetic field strength (panel (a)), temper- 
ature at logr = 0.5 and logr = —2.0 (panels (b) and (c)), and 
temperature gradient offset Torad as defined in Equation [y (panel 
(d)). The contour lines of the magnetic field map (panel (a)) show 
the temperature at logr = -1.0 (5200, 5500, and 5800 K), the con- 
tour lines in the lower three panels show the height-independent 
magnetic field strength (250, 500, and 1000 G). 

magnetic patches. 

Measured (symbols) and fitted (lines) Stokes I and V 
profiles within a network patch (pixel 1, red), a gran- 
ule (pixel 2, blue), and an intergranular lane (pixel 3, 
yellow) are shown in Figure |4] (pixel numbers (764,809), 
(667, 805), and (713, 742) in the IMaX data set 163-209). 
Owing to its low excitation potential of 0.12 eV, the Fel 
5250.2 A line is highly sensitive to temperature. The high 
temperature at the line formation height in the magnetic 
patch (1) therefore results in an extremely weak Stokes / 
profile, with a line depth of only 6% of the local contin- 
uum (see Figure HJ. At the same time, the continuum is 
higher than the average quiet- Sun continuum. The line 
depths of a typical granular profile (2) and an intergran- 
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Fig. 4. — Measured (symbols) and fitted (lines) Stokes I (top) 
and V (bottom) profiles for the pixels (1), (2), and (3) in Fig- 
ure [T] The vertical, dotted lines indicate the IMaX filter positions. 
The decrease in continuum intensity at both boundaries and the 
increase in Stokes V at the red boundary are caused by the neigh- 
boring Co I and Fe I lines. 
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Fig. 5. — Temperature profile for model atmospheres and as re- 
trieved from the inversions for the pixels (1), (2), and (3) in Fig- 
ure [T] 

ular profile (3) are significantly larger (~57% and 56%, 
respectively, of their continuum levels). 

In Figure [5j we compare the temperature stratifications 
returned for the three selected pixels featured in Figure 2J 
The temperature gradient in the intergranular lane (3) is 
somewhat flatter than of the HSRASP atmosphere. The 
temperature in the center of the granule (2) is higher 



in deep layers and falls off rapidly with height. The 
kilo- Gauss network patch (1) shows high temperatures 
that, for the m ost part, lie between the empirical network 
JSoLanki 19861. labeled NET (S86)) and the plage models 
(jSolanki fc Bridieviclll992L PLA(SB92)). In the upper 
layers it is even hotter than both of these atmospheres, 
but not well constrained. Similar temperature strati- 
fications (and also field strengths) were obtained from 
the analysis of spectropolari metric observations in the 
Fel 6 30 nm line pair (see, e.g., iSanchez Almeida fc Litesl 
l2000h . but now, for the first time, we obtain this tem- 
perature stratification without introducing a magnetic 
filling factor. The relevant atmospheric parameters re- 
sulting from the inversion for pixel (1) are To = 6467 K, 
T G rad = -652K/logr, B = 1447 G, 7 = 169°, and 
£ mic = O.lkms- 1 , for pixel (2) T = 6865K, T Grad = 
167K/logr, and £ m i c = 1.4kms _1 , and for pixel (3) To = 
6139K, T Gr ad=-168K/logr, and £ mic = l^kms" 1 . 
The inclination of the magnetic field for pixel (1) is likely 
to be more vertical to the solar surface since a small linear 
polarization signal could not be well constrained by the 
inversion. The kilo- Gauss nature of the magnetic patch 
persists even under the assumption of a completely verti- 
cal (180°) field. The magnetic field strength and inclina- 
tion for pixels (2) and (3) cannot be determined reliably 
due to the low polarization signal. 

5. SUMMARY AND CONCLUSIONS 

The IMaX instrument on board the Sunrise mission 
allowed photospheric magnetic fields to be analyzed with 
unprecedented resolution. The quiet-Sun region stud- 
ied in this Letter contains two network patches with 
the characteristic observational signature of an extremely 
weak Stokes / profile, associated with a relatively strong 
signal in Stokes V. 

Without requiring a magnetic filling factor, the 
inversion of the Stokes profiles of a single pixel 
within such a network patch returns a magnetic field 
strength of 1.45 kG. This value is consistent with field 
strengths determined from the Fel 5250.2 A line in the 
quiet-Sun network usin g indirect techniques. Thus, 
IStenflo fc HaTv^vl (Il985h find va l ues o f 800-1100 G. Sim- 
ilarly, iGrossman n-Doerth et al.l (fl99 6) obtained a distri- 
bution of intrinsic field strengths in the quiet-Sun from 
this line, which peak around 800-1000 G. 

The extremely weak Stokes / profiles in the network 
patches point to their high temperature in the higher at- 
mospheric layers. We find the temperature at heights 
of logr = —1, corresponding to the formation height 
of the line core in these areas, to lie significantly above 
the temperature in the surrounding area. The inver- 
sions directly give rise to an atmospheric stratifica- 
tion of physical quantities consistent wit h existing flux 
tube models in network plage regions (jSolankil Tl986t 
iSolanki fc Brigljevid Il992f ) , and with the properties of 
magnetic flu x concentrations i n magnetohydrodynamic 
simulations (|Vogler et al.ll2005[ ). 

We can identify 5-10 similar network patches within 
every IMaX snapshot. Inversions of these network 
patches yield field strengths of up to 1.8 kG. The tem- 
perature stratification and the other atmospheric param- 
eters are similar to the values presented in Figure [5] and 
in Section [4] for pixel (1). Different from the analysis 
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of iBerger et al.l ([2004f h these kilo- Gauss structures are 
no obvious flux concentrations seen as micropores in the 
intensity images. 

Both results, the kilo-Gauss field strength and the high 
temperatures in the mid- to upper photosphere, lead us 
to the conclusion that the observed network patch is a 
magnetic flux tut>4E The unprecedented spatial resolu- 
tion of achieved with IMaX on Sunrise, was 
sufficient to resolve this flux tube. This has been a long 
cherished aim in solar physics. We have demonstrated 
that with aim class telescope under ideal conditions it is 
achievable, even in the quiet-Sun. The data set obtained 
during the Sunrise flight will allow us to analyze the 
temporal evolution of these flux tubes in the quiet-Sun 
as well as to study their properties in a statistical sense 
without dilution due to serious intermingling with light 
from field- free gas. Additionally, the connection between 
the shallow temperature gradient in these flux tubes and 
the bright points at chromospheric heights observed with 
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1 The rapid expansion of a vertical flux tube with height re- 
sults in a strong magnetic field gradient, which was not included 
in this analysi s. Test inversions including a magnetic field gradient 
of -3 G km -1 (jMartmez Pillet et "811119971 ) reproduce the observed 



Stokes profiles equally well as inversions with a height-independent 
magnetic field strength. This gradient model yields field strengths 
at the formation height of the line core (log r ~ 1) very similar to 
the field strengths presented in this Letter. 



